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Abstract 

This letter reports the first detailed measurement of M z , the component parallel to the effective 
field direction, when ferromagnets are excited by microwave fields at high power levels. It is found 
that M z drops dramatically at the saturation of the main resonance. Simultaneous measurements 
of M z and absorption power show that this drop corresponds to a diminution of the spin-lattice 
relaxation rate. These changes are interpreted as reflecting the properties of longitudinal spinwaves 
excited above Suhl's instability. Similar behavior should be expected in spinwave emission by 
currents. 
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The out of equilibrium dynamics of ferromagnets is receiving much attention owing to the 
potential application to spin electronic devices. Interest for transport properties usually lies 
in the response of the longitudinal magnetization M z (component parallel to the effective 
field direction) to a large excitation amplitude. But the high power dynamics cannot be 
understood in terms of a linear susceptibility. Nonlinear (NL) contributions contained in 
the torque term of the gyroscopic equation become important as soon as the precession angle 
exceeds a couple of degrees. Although these effects were first discovered in the ferromagnetic 
resonance (FMR) response of insulators , they are in fact generic to all ferromagnets. The 
same physics also applies to metallic samples as shown recently by An et al. Qj. While the 
consequences of these nonlinearities on the microwave susceptibility have been thoroughly 
investigated j^, their effects on M z have never been established. 

This letter reports the first measurement of M z at the resonance saturation. We will 
present our data obtained at room temperature on an yttrium iron garnet (YIG) sample, 
shaped into a disk of diameter D = 160/im and thickness 4.75/im and uniformly magnetized 
by a static magnetic induction B ext applied along the normal axis z of the disk 4| . It is found 
that M z drops dramatically at the saturation of the main resonance. It results from a rapid 

[rowth of phase-coherent spinwaves (SW) propagating along the magnetization direction 
|. We use these experimental results to discuss the consequences of the NL effects on 
theproblematic of spin injection. Their inclusion inside the microscopic spin-torque model 
IfL 1 should provide a mechanism to reverse the magnetization through a reduction of 

m z y,y. 

We exploit the exquisite sensitivity of Magnetic Resonance Force Microscopy (MRFM) 
to follow the changes in M z with a resolution better than lppm. A schematic of the setup 
is shown in Fig^a). A cylindrical permanent magnet is glued at the free end of a clamped 
cantilever and then aligned with the axis of the YIG disk. The distance between the sample 
and the probe is fixed at 100/im so that their coupling is in the weak interaction regime 
The inhomogeneous dipolar field generated by the sample creates a point load on the tip 
and thus an elastic deformation of the cantilever. Since the flexural modes of the cantilever 
are well below the Larmor frequency, the mechanical probe is insensitive to the precession 
of the transverse magnetization and it only couples to the longitudinal component M z . 

The measured quantity is AM 2 = M s — M z when the microwave field is turned on. M s 
is the saturation magnetization at the lattice temperature. The microwave is generated 
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by a synthesizer and fed into an impedance matched stripline resonator tuned at ujq/2-k = 
10.47GHz. We call h the circularly polarized amplitude of the microwave driving field 
at the sample position. If h is much lower than Suhl's threshold (see below) then the 
response is proportional to the excitation power. Fig^b) shows the FMR spectrum of our 
disk in the linear regime. A multiplicity of AM Z maxima are detected during a sweep of 
B ext . The spectrum illustrates the quantized energy levels of magnetostatic waves (near- zero 
wavevector) confined by the sample diameter. 

Simultaneously, we measure the imaginary part of the transverse susceptibility x" through 
a standard setup. The power reflected off the half-wavelength resonator is detected by a 
microwave crystal diode whose signal is proportional to the absorbed power P a b s = wox"h 2 . 

Although the discussion below can be extended to any of the magnetostatic modes above, 
we concentrate on the high power behavior of the main resonance at B = 5324. 5G (see 
Fig^b)). This mode has a transverse wavevector k pa ir/D and is called hereof the uni- 
form precession because it has no precessional nodes. We plot in Fig j2] the microwave field 
dependence of x" (closed circles) and AM z /h 2 (open circles) evaluated at resonance. Both 
quantities are normalized by their low power values. The x" data shows the well known 
premature saturation behavior above Suhl's threshold, h c = 5mOe. Surprisingly, the power 
dependence of AM z /h 2 exhibits a peak at the threshold. Evidence of a behavior where 
AM z /h 2 increases with power has never been established before in any magnetic resonance 
experiment. 

Obviously the correlation between the longitudinal and transverse quantities cannot be 
explained by a macrospin picture, which implies that AM Z = x" 2 h 2 /(2M s ) in this power 
range. The difference clearly demonstrates that the x" measurements only reveal a partial 
picture of the dynamics inside the sample. Other SW with identical energies (ou^ = Uq) 
but different wavevectors (k ^> k ) can be excited in the sample. FigJ3fa) is the magnon- 
manifold for our finite aspect ratio disk [ijj. The shaded area expresses the possible balance 
between demagnetizing and exchange energies. Volume demagnetizing energy of SW mode 
can be changed by varying the angle 8k between k, the propagation wavevector and M s , 
the magnetization vector. Short wavelength SW do not couple to x" because the transverse 
projection averages to zero. However, each magnon excitation diminishes the longitudinal 
component by 7ft. In other words, while the transverse data measures the number of uni- 
form magnons n u = \ X M h h , the longitudinal data is proportional to the total number of 
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magnons n t = Here we only count the magnons created by the microwave excitation. 

Thermal magnons' occupation numbers are accounted for in M s and remain constant in all 
the measurements below. Suhl showed jj| that the NL coupling destroy the independence 
between degenerate modes. He found that, above a critical power (h^), n u saturates at n c 
as shown in FigEfb). We want to establish that this saturation is associated with a surge 
of degenerate magnons nt — n u . 

Before continuing further the discussion, we need to look at the different relaxation chan- 
nels in our sample. These damping forces are necessary to stabilize the SW excited by the 
nonlinearities |3j . Fig0] establishes a schematic of the coupling between the various degrees 
of freedom in the system. At ujq and Bo, the uniform microwave field preferentially couples 
to the uniform mode, the longest wavelength mode available inside the degenerate band. Its 
energy relaxation rate directly to the lattice will be written 770 (our definition is twice as 
large as the amplitude relaxation rate used in Suhl's paper [3]). The linewidth of the peaks 
in FigGJb), however, is the sum (t/ + ^sp)/7, where r) sp is the decay constant of the uniform 
precession to degenerate SW due to scattering on the sample inhomogeneities. We define 
as the number of degenerate magnons having a wavevector k and 1% as their average decay 
rate to the thermodynamic equilibrium. Translated into Bloch's notation, the linear part 
of the diagram in FigJU corresponds [3| to T 2 = 2/(770 + r] sp ) and T\ = T 2 /2 (1 + rj sp /rjk), 
respectively the transverse and longitudinal relaxation times. These relaxation times, in the 

imit of infinitesimal excitations, have been completely characterized in a previous paper 

1^ | . Their values are summarized in Table H] 

The NL terms in the gyroscopic equation Q] couple coherently degenerate SW of equal 
and opposite wavevectors (+k, —k). The coupling term with the uniform motion, £fcn u , is 
maximum (£fc| max /7 = 2nM s ps 900G) for longitudina SW (8 k = or fc||M s ). We use 
Suhl's notation for so that the reader can refer to j^l for a complete treatment. In our 
disk, these SW have a wavevector fc max 6.3 x 10 4 cm _1 (see Fig|3fa)). A critical threshold 
is reached when £fcn u becomes comparable to r] z , the relaxation rate to the lattice of the 
longitudinal SW. Two quanta of the uniform precession then break down in a parametric 
magnon pair propagating along M s . The instability corresponds to a spatial distortion of 
the instantaneous axis of precession which diminishes the transverse demagnetizing energy. 
The critical number of uniform magnons is n c = | lMa J 2 hc with h 2 c = ^2^7^, the saturation 
power. We propose to view this as an alteration of the energy flow between the modes in 
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Figg] 

w;n u = - {t] + r] sp f(n u )}n u (la) 
d 

j^n x = r] sp {f(n u ) - l}n u - r] z n z , (lc) 

where f(n u ) = 1/ \/l — n^/n 2 jlJ] and n z is the number of parametric magnons. 

In the stationary regime, EqlTal leads to an implicit equation for the susceptibility first 
derived by Suhl |l4l |: 

? = "; + "i 4 (2) 

Vo + Vsp/yl-x" {h/h c f 
where x" = x" I ilT 2 M s ) is the susceptibility normalized by its low power value. We try 
to compare quantitatively our data to the profile predicted by EqEl We use the measured 
value of the threshold, h c = 5mOe, in the formula. It corresponds to 772/7 = 0.15G, a value 
which is consistent with previous measurements in spheres and t hin p lates where it was 
found that the decay rate becomes smaller as 9k diminishes [12I We have displayed in 
Fig|2]the predicted behavior, with no fitting parameters. The model properly predicts the 
shape and the smearing out of the singularity at h c . 

Solving the set of Eq^in the steady-state gives a new analytical relationship between x" 
and AM z /h 2 

AM, 1 r~ ~2i ~2 



v 21,,, i*" -*")+*"■ < 3 > 

where AM Z = AM z / \<y 2 TiT 2 M s ) is normalized by the slope of AM z /h 2 at low power. The 
above equality is equivalent to an energy balance between the absorbed power and the 
dissipated power in the sample P a b s = -Pdiss = ^0X4/0}%^ Q, Qj- The solid line in 
Figgis the bell-shape behavior inferred from EqOlwhen Tirj z = 0.15. Experimentally we 
find that the raise of AM z /h 2 begins at powers lower than h c . We have omitted in our 
model the finite size effects in order to keep the discussion simple. But the so-called uniform 
mode has a spatially dependent precession amplitude. AM Z at the center of the disk is 



rage J|. Changes in the profile are expected 
when different locations of the sample hit the saturation threshold. Spatial dependences are 
further magnified by the MRFM which provides a local measurement (here the mechanical 
probe is above the disk center). While the raise of AM z /h 2 coincides with the saturation 



of the central part of the disk, the spatially averaged transverse susceptibility x" drops at 
higher power, approximately when the periphery saturates (zone with the highest spatial 
weight). The dashed line in FigElis the predicted behavior when EqEJis weighted by the 
spatial profile of AM Z . Such good agreement with the data suggests that the shift of the 
peak towards h < h c is due to these finite size effects. 

To gain further insight, we plot in Fig|S]the ratio between the transverse and longitudinal 
quantities. This shows the power dependence of the spin-lattice relaxation rate inside the 
spin system. In magnetic resonance, the power absorbed during the spin-lattice relaxation 
time is equal to the energy stored in the sample, ^AM Z . Thus the quantity on the ordinate 
axis is ^j^f- = ^ 12{k} Vk^- We find that the observed correlation between x" an d AM z /h 2 
corresponds to a monotonic diminution of the damping with increasing power. This ratio 
starts to decrease below h c and the drop extends to the h > h c region. In our model, the 
relaxation rate is simply equal to 1 /T\ when f(n u ) = 1 (linear regime) and then it reduces to 
rj z when f(n u ) ^> 1. The two limits correlate respectively to n t ~ n u and n t ~ n z . The solid 
line in FigEl is the drop inferred from the analytical expression showing how the r] z value is 
approached asymptotically. We find that the agreement with the data extends well above 
h c and the asymptotic value is in agreement with the r\ z inferred from the threshold value, 
h c . We recall that our model only considers the NL coupling between the uniform and the 
longitudinal magnons. A proper analysis, however, should also account for the NL coupling 
between all the degenerate modes, which become important at much higher powers. They 
correspond to a different redistribution of the degenerate magnons' occupation number in 
the flow diagram of FigEJ 

In conclusion, this letter studies the changes of properties of excited SW due to NL 
effects. If the excitation mechanism favors a particular mode, then a 4-magnon process starts 
at high power levels to redistribute the magnons' occupation number between degenerate 
modes. During the redistribution, the coupling to the excitation usually decreases: less 
energy is absorbed for the same excitation energy. The efficiency, however, might increase 
if the newly filled modes (generally the longitudinal ones 111 ]) have lower energy relaxation 
rates: more SW are emitted for the same absorbed energy. We believe that these comments 
are also pertinent to SW emitted by a current for the following reasons. First, the above 
statements can be generalized to other selection rules, i.e. the favored mode might not 
necessarily be the longest wavelength one. Second, the precession amplitudes achieved by 
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currents are much larger than the ones achieved by microwaves |l7j (though conductors 
have larger values of 1] ). Third, the eigen-modes are at least doubly degenerate if one 
sample dimension is larger than lOOnm (indeed the width of the magnon manifold is of the 
order of 27r//c max < lOOnm for metals). Fourth, the energy relaxation rate of degenerate 
SW generally decreases as k aligns with the M s direction (cf. measurements in Permalloy 
films j^]). These arguments indicate that parametric SW should also get excited by spin 
injection. Excitations of short wavelength SW should provide a mechanism to reverse the 
magnetization through a reduction of M z Finally, we emphasize that the NL transfers 
conserve the phase-coherency. This point is important regarding to the spin torque model 
which relies on the concept of coherent emission of SW p] . 
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FIG. 1: (a) Schematic of the experimental setup. A Magnetic Resonance Force Microscope is used 
to measure AM Z while a standard microwave setup measures x" ■ (b) Linear FMR spectrum of the 
micron-size YIG disk detected mechanically. The main resonance occurs at Bq. 
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FIG. 2: Microwave field strength dependence of the transverse (closed circles) and longitudinal 
(open circles) components of the magnetization at Bq . The solid lines are the analytical predictions 
(see text). The dashed line is the behavior corrected by finite size effects. 




FIG. 3: (a) SW dispersion: the shaded area is the magnon manifold explained in ll| (A is the 
exchange constant), (b) Power dependence of n u and nt — n u for our sample. 




excitation 



uniform degenerate 
precession (ko ^ 0) niagnons (k ^ 0) 

n u n t — n u — n z + n. 



lattice 



phonons 



thermal 
magnons 



FIG. 4: Block diagram of the energy transfers between the various degrees of freedom of the system. 
Transfers between the upper reservoirs conserve the total energy. The linear regime corresponds 
to f(n u ) = 1. The shaded area and doubled line are the alterations by NL effects. 
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FIG. 5: Power dependence of the ratio of x" over AM z //i 2 . It measures the energy decay rate of 
the spinwave system. The solid line is the analytical prediction. 
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mh 
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Vz/l 


0.2 


1.07 


0.65 


0.15 



TABLE I: Relaxation rates of degenerate magnons (in G). 
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